A Quintom scenario with mixed kinetic terms 



Emmanuel N. Saridakis 1 E and Joel M. Weller 2 'Q 

1 College of Mathematics and Physics, Chongqing University of Posts and Telecommunications, 

Chongqing, 400065, P.R. China 
2 Department of Applied Mathematics, University of Sheffield, 
Hounsfield Road, Sheffield S3 7RH, United Kingdom 

We examine an extension of the quintom scenario of dark energy, in which a canonical scalar field 
and a phantom field are coupled through a kinetic interaction. We perform a phase space analysis 
and show that the kinetic coupling gives rise to novel cosmological behaviour. In particular, we 
obtain both quintessence-like and phantom-like late-time solutions, as well as solutions that cross 
the phantom divide during the evolution of the Universe. 
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I. INTRODUCTION 



According to independent observational data, includ- 
ing measurements of type la supernovae the Wilkin- 
son Microwave Anisotropy Probe @], the Sloan Digital 
Sky Survey [|| , and X-ray observations Q , the Universe 
is undergoing a phase of accelerated expansion. Although 
the cosmological constant seems to be a simple and eco- 
nomic way to explain this behaviour [BJ , the extreme de- 
gree of fine-tuning associated with this means that it is 
attractive to consider other dark energy candidates. The 
dynamical nature of dark energy, at least at an effective 
level, can originate from various fields, such as a canon- 
ical scalar field (quintessence) @ or a phantom field Q 
(i.e. a scalar field with negative kinetic terms). Another 
possibility of considerable interest is the quintom scenario 
[8l4l0(. in which the acceleration is driven by a combina- 
tion of quintessence and phantom fields. The quintom 
paradigm has the advantage of allowing the dark energy 
equation-of-state (EOS) parameter (u>de) to cross the 
phantom divide during the evolution of the Universe, an 
intriguing possibility for which there is observational ev- 
idence from a variety of sources (see [13 for a review). 

The aforementioned models offer a satisfactory descrip- 
tion of the behaviour of dark energy and its observable 
features. However, the dynamical nature of dark energy 
in these scenarios leads to the "coincidence" problem, 
which many authors have sought to resolve by consider- 
ing a coupling between dark energy and the other com- 
ponents of the Universe. Thus, various forms of interact- 
ing dark energy models have been constructed, including 
coupled quintessence (Tl| and interacting phantom mod- 
els [121 ]. Similarly, the dynamical nature of dark energy 
makes it difficult to fulfil the observational requirement, 
w de ~ — 1 [HI, which has led to the introduction of ad- 
ditional scalar fields. Inspired by similar multi-field infla- 
tion [l4| and assisted inflation constructions [l5[ , various 



forms of assisted dark energy models have been consid- 
ered [IE Ell- 
One generalised class of the multi-field quintessence 
scenario allows for a mixing of the fields' kinetic terms 
[lH (see also [lj|). Mixed kinetic terms appear also in 
the non-trivial scalar field models based on string theory 
constructions, studied in the framework of generalised 
multi-field inflation [2fj| |. 

In the present work we aim to combine the advan- 
tages of the quintom scenario with those of assisted 
quintessence with mixed kinetic terms. Thus, we con- 
struct a quintom model in which the kinetic terms of the 
canonical and the phantom fields are mixed. Indeed, as 
we see, the model exhibits features characteristic of both 
quintessence and phantom models. 

The plan of the work is as follows. In Sec. |TT] we 
construct the quintom cosmological scenario with mixed 
kinetic terms and present the formalism for its transfor- 
mation into an autonomous dynamical system. In Sec. 
IIIII wc perform the phase-space stability analysis and in 
Sec. IIVI we discuss the cosmological implications of the 
results. Finally, our conclusions are presented in Sec. [V] 



II. QUINTOM SCENARIO WITH MIXED 
KINETIC TERMS 

Let us construct the quintom cosmological scenario 
with mixed kinetic terms. Throughout the work we con- 
sider a flat Robertson- Walker metric: 



ds 2 



dt l 



a 2 (£)dx 2 



(1) 



with a the scale factor. 

We consider a model consisting of a canonical scalar 
field 0, a phantom field a and a mixed kinetic term pro- 
portional to the parameter a. Thus, the action, in units 
where 8irG = 1 reads: 



S = 
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g^d^d v a + V^) + V a {a) + S M , (2) 
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with V^(4>) and V a (cr) the corresponding field poten- 
tials, and Sm the action for the matter (dark plus bary- 
onic) component of the Universe. We assume that the 
matter is described by a perfect fluid with energy den- 
sity pm, pressure pu and barotropic index 7. Thus, 
Pm = (7 — l)/"Af an d the matter equation-of-state pa- 
rameter is wm =7—1. Varying the action with respect 
to the metric yields the Friedmann equations 



H 2 = 
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T 2 
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+V (j> {4>) + V (T {a)+ PM 
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-V^) - V a (a) +p M 



(3) 



(4) 



while the evolution equations for the two fields read 



da 2 d<f> ) 



"^T + i^ U + ^T =0- (6) 
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In quintom scenarios the dark energy is attributed to 
the combination of the two scalar fields. In the case at 
hand, from the Friedmann equations (J3J) , (J4j we straight- 
forwardly read the dark energy density and pressure as 

Pde = \^-\& 2 + ^& + V^)+V a (a) (7) 
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Pde = ^--&' + -ct>a-V^)-V a (a), (8) 

and thus for the dark energy equation-of-state parameter 
we obtain 



Pde = \<\> 2 - \° 2 + § <l>& ~ W) ~ V *(v) 
Pde 



102 _ 1 & 2 + + V ^ + 



(9) 

Finally, using jZ]),® and ([5]), ([6]), we can easily acquire 
the dark energy evolution equation in fluid terms, which 
takes the expected form 



Pde + 3H(p DE + p DE ) = 0. 



(10) 



We mention that in the case a = the above scenario 
coincides with the usual quintom one @. 



III. PHASE-SPACE ANALYSIS 

In order to perform the phase-space and stability anal- 
ysis of the quintom model at hand, we have to trans- 
form the aforementioned dynamical system into its au- 
tonomous form (2l| - [23j . This can be achieved by intro- 



ducing the auxiliary variables: 



£0 



s/6H 



together with N = In a. It is easy to see that for every 
quantity F we acquire F = H^. Using these variables, 
from ([7]) we can calculate the dark energy density param- 
eter, 



n _ Pde 2 2 , ,2,2 



t-DE 



3H 



2 X (p X a 



ax^Xa + y <t , + y a , (12) 



and from ([9]) the dark energy equation-of-state parame- 
ter, 



xl-xl + ax^x a -yl-yl 

W DE = 3 _ n ,2 , I ..2 , ..2 ' ( 13 ) 



The final step is the consideration of specific poten- 
tial forms. Following works on assisted quintessence fl7T | 
and assisted inflation 15| scenarios, to geth er with vari- 
ous other cosmological models [U |24| , we consider 
exponential potentials of the form, 



(14) 



Note that equivalently, but more generally, we could con 
sider potentials satisfying A — — dv W 

1 dv( ■ 



const, 



KV(tp) d</> 

^ = ~ kVTo) ' da ' ~ cons t (f° r example this relation is 
valid for arbitrary but nearly flat potentials [25|). Fi- 
nally, without loss of generality, in this work we consider 
A and p to be positive, since the negative case alters only 
the direction towards which the fields roll, and can al- 
ways be transformed into the positive case by changing 
the signs of the fields themselves. 

Let us now perform the phase-space analysis of the 
model. In general, having transformed the cosmological 
system into its autonomous form X' = f(X), where X 
is the column vector constituted by the auxiliary vari- 
ables, f(X) the corresponding column vector of the au- 
tonomous equations, and prime denotes a derivative with 
respect to N = In a, the critical points X c are extracted 
satisfying X' = 0. In order to determine the stability 
properties of these critical points we expand around X c , 
setting X = X c + U with U the perturbations of the 
variables considered as a column vector. Up to first or- 
der we acquire U' = Q • U, where the matrix Q contains 
the coefficients of the perturbation equations. Thus, for 
each critical point, the eigenvalues of Q determine its 
type and stability. 

Following this method and using the auxiliary variables 
(fTTj) . the cosmological equations of motion Q, @ 
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and (|6|) become 
x'j, = -3au + 



-3x a 



3 f a 2 

2V 1 + T 

3 / a 2 

2V 1 + T 



Ay, 



Wo 



Va = 

where 



-Xx^y^ + y$T, 



-^x a y a + y a T, 



Hy~ j + x^T, 
(15) 



2x 2 - 2x1 + 2ax^x 



+7(1 ~x\ + x\- a X(t> x a -y\- yl) . (16) 

The critical points (x^ c ^x ac ,y^ c ,y a ^) of the au- 
tonomous system (fT5]) are obtained by setting the left 
hand sides of the equations to zero. The real and phys- 
ically meaningful (i.e. corresponding to y^^a > 0, and 
< Qde < 1) points are presented in Table |TJ 

The 4x4 matrix Q of the linearised perturbation equa- 
tions of the system (|15[) is shown in the Appendix. Thus, 
for each critical point of Table [J we examine the sign of 
the real part of the eigenvalues of Q to determine the 
type and stability of the point. In Table |TT] we present 
the results of the stability analysis. In addition, for each 
critical point we calculate the values of SIde [given by 
(IT2"|) ] and of wde [given by ([15])]. Finally, we mention 
that the labelling of the critical points follows the con- 
vention used in [18l |. That is, since in a two-field system 
there are critical points in which the potential term of one 
of the fields is zero, the connection between these points 
is made more explicit. For instance, in the points Dl and 
D2 (for which Q,de = 1) only one of the fields' potential 
terms contributes to the energy density. Note that in this 
model (unlike that in pj|) the solution satisfying X' = 
with j/0 = corresponding to CI is unphysical. 

The stability of the critical points can be summarised 
as follows: 

• Point A 

This is a trivial solution,corrcsponding to the fluid 
dominated point F in |8|, where the kinetic and 
potential components of both the phantom and the 
quintessence fields are negligible. It exists for all 
values of a, A and /1. The eigenvalues are 



ei 



§7>0, 



e2 

e 4 = -f (2- 7 ) <0, 
so this is a saddle point. 
• Point B 



(17) 



The two points indicated by the ± are similar to 
point B in [Hi. When a = (point Kin 0), there 



is a continuous locus of points that describes a hy- 
perbola in the {x^.Xg-) plane. When a^O the lo- 
cus still describes a hyperbola, but the asymptotes 
are no longer orthogonal. The eigenvalues are 

3(2- 7 )>0, 

0, 

q _ 5^6 , 



ei 

f'3 



3V6 
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1 ^ 



:± \yjxl((x 2 +4) +4 



4 p x o 
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\x a a±\ V4(« 2 +4)+4j , (18) 

where the ± in e% and corresponds to the par- 
ticular solution considered. e\ is positive, therefore 
the point is unstable. 

• Point Dl 

This takes a similar form to point Dl in (The 
value in the solution corresponding to CI in [l8T | 
has a factor of \/(4 + a 2 )(l — 2/7) so there are no 
real values of a for which the point exists.) In the 
standard quintom case (a = 0), Dl is a phantom- 
dominated point (point P in the notation of 0]). 
In that case, as in the present one, the term under 
the square root in y ac is positive and Qde = 1. 
so the point exists for all values of the parameters. 
However, in the case at hand, the eigenvalues read 



ei,2 



B3 



(4 + a 2 

4/z 2 
(4 + a 2 
fi(2fi 



3 < 0, 



3 7 <o, 



aX) 



(4 + a 1 ) 

so the point is stable only if (2/i + aA) > 0, which 
is trivially satisfied when a > 0. 

• Point C2 

The points C2 and D2 are similar to Dl, in that 
only one of the fields' potential energy contributes 
to the total energy density; in Dl this was the 
phantom field, while here it is the quintessence 
one. This means that there is both a scaling so- 
lution (C2) and a field-dominated solution (D2), 
for which Ude = 1- The density parameter for C2 



is n DE = 4^(4 

means that A > 



a ), so the condition Ude < 1 



f (4- 



1 a 2 ) must be satisfied if 
the point is to be physical. The eigenvalues are 



ei 

('2 

X 



-|(2-7)<0, 



Q/x), 



3(2-7) 
4 



1± 



'6 7 2 (4 + a 2 ) + A 2 (2-9 7 ) 
A 2 (2- 7 ) 



(20) 
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V<t>c 
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A 
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-±x a a± iy'iJ(o ;i +4)+4 


3-D" 








Dl 













3(4+0^) 


3(4+a 2 ) 


3(4+a^) 


C2 


iVfi 
2A 


Q7V6 
4A 


^(4 + ^)(|_ 1) 





D2 


2^6 A 
3(4+o^) 


v'BAa, 
3(4+o^) 


3 P ^ 7 ^(4 + a2 )(6+ 3 Q 2_ A 2) 





E 


V 2 A 


V 2 K 




Vl^^(fM-A)(f-l) 


F 


n 

76A 


n 


TTA^^+fA) (l-f) 




nA(f M - a) 




TABLE I: The real and physically meaningful critical points of the autonomous system (|15[). The parameter II is defined in 
J22J. 




Existence 


Conditions for Stability 




WDE 


A 


all fi,\,a 


unstable 





undefined 


B 


all fi,\,a 


unstable 


1 


1 


Dl 


all /i,A,a 


fj, > —aX/2 


1 


1 4 M ^ 

3(4+a a ) 


C2 


A>^(4 + o») 


A < a/i/2 


4^(4 + « 2 ) 


WM 


D2 










1 


1 1 4a2 


A< V / |(4 + a2 ) 


\<min an/2,sj^(4 + a 2 ) 


1 ^ 3(4+cr*) 


E 


QAt - 2A > ; 2/i + qA > ; II > 37 


unstable 


37/n 




F 


a/i - 2A > ;0<n<6 
+ <0 


unstable 
stable 


1 


-i + n/3 



TABLE II: The properties of the critical points of the autonomous system (|15[) . The parameter II is defined in (122[) . Note that 
there are two disconnected regions in the parameter space for which point F is physical. 



e-i and are always negative as the condition for 
existence gives 

6 7 2 (4 + a 2 ) + A 2 (2 - 9 7 ) < A 2 (2 - 7). 

Thus, the point is stable if (2A — a/i) < 0, so un- 
like the quintom case (cf. point T in Q), with a 
sufficiently large value of afi one obtains a stable 
scaling solution in this system. 

• Point D2 

This field-dominated solution is similar to Dl in 
[HI and reduces to point S in Q when a = 0. The 

existence condition, A < y^|(4 + a 2 ), is derived 

from the requirement that y$ c € R. The eigenval- 
ues are 



ei 



A(2A - an) 
(4 + a 2 ) 



4A 2 



''2 



2A 2 



e3,4 



(4 + a- 2 



(4- 



-3 < 0, 



3 7 



(21) 



where the existence condition ensures that 63.4 are 
negative. If the point is to be stable both e\ and 



e.2 must be negative, which is satisfied if 



A < min 



J ^ (4 + a 2 ) 



Once again, one can choose values of \i and a so 
that this point is stable, which is not the case in 
the standard quintom model. 

• Point E 

The points E and F are characterised by the fact 
that the potential terms for both the phantom and 
the quintessence field are non-zero. We mention 
that when a — these points do not exist, since 
it is not possible to choose values for A and n such 
that y^cVac € R- Point E is a scaling solution, 
with density parameter 



n 



DE 



37 

AV 



(/i 2 — A 2 + a/xA) 



Following [18j , we can rewrite this quantity in terms 
of an effective exponent II (which in the present 
case can be negative too, contrary to [l~8[) defined 



5 



by 



The eigenvalues read 



1 

IT 



1 



I'- 



cx 

jj,X ' 



(22) 



which leads to the existence condition II > 37. 
Note that, as /i, A > 0, 



{afi - 2A) > ^> (2/i + aX) > 0. 



(23) 



Thus, to ensure that y^ c and y ac are real, one must 
additionally impose (a/i — 2A) > 0. The eigenvalues 
are 



ei,2 



e3,4 



1(2-7) 



1 ± 



'247 2 /Tl + 2-97 



(2-7) 



■1(2-7) 



1± 4 1 



87(2/1 + aA)(2A-a/i) 
A/i(2- 7 )(4 + a 2 ) 



(24) 



(25) 



In terms of II the existence condition reads 



247^ 



2 - 97 < (2 - 7); substituting this into (HH leads 
to ei, 2 < — 1(2 — 7)(1 ± 1) < 0. However, looking 
at the fraction under the square root in ([2"5]), we 
observe that 



87(2/1 + a\){2\ - afi) 
A/i(2- 7 )(4 + a 2 ) 



> 1. 



Thus, either e 3 or must be positive, and therefore 
the critical point is unstable. 

• Point F 

In this (field-dominated) point we observe the pres- 
ence of square-root terms in y$ c and y ac . Thus, for 
these values to be real we require 



(2/z + aA) > 0, 

(a/x-2A) > 0. 

Using ([22]) and (f23]> . and noting that 
(2/i + aA) < (a/i - 2A) < 0, 

and 



(2/x + aA) < 



fi — A" + a/(A < 
n < 0, 



(26) 



(27) 



we conclude that there exist two possible cases: 

(a) a/i-2A>0; < II < 6, (28) 
(6) 2/i + aA < 0. (29) 



-3 1 



ei 



e 2 = n - 37, 



3 / n 

e 3A = -2(1-6 



1± 4/1 



4n (2/i + «A)(2A-a/x) 
3V (4 + a 2 )(l-S) 



(30) 



In case (a), (l — §) > 0, thus e\ < 0. However, in 
a similar way to the analysis of point E, the term 
under the square root in the expression for 63.4 is 
greater than 1 so either £3 or is positive and the 
point is unstable. In case (b), the term under the 
square root is less than 1 so 63,4 < 0. ei < is triv- 
ially satisfied [using (|2"7j)] and, again, (l — S) > 0, 
therefore e\ < 0. Thus, the point is stable whenever 
the condition (|29p is satisfied. Finally, we mention 
that this separation of the parameter space is a 
novel feature of the scenario at hand, which is not 
present in the corresponding case of two canonical 
fields 

The stability regions for the points are plotted in Fig. 
[TJ In the standard quintom case (a = 0), there is only one 
stable critical point that exists for all values of \i and A, 
corresponding to the phantom point Dl. In the case a < 
the parameter space is fractured, thus, although there 
is a stable point with wde < — 1 at each parameter-space 
region, the precise value of »de for particular values of 
fi and A depends on whether point Dl or F is stable. In 
the case where a > 0, the phantom dominated point Dl 
is always stable. It is interesting, however, that there 
are also regions where C2 (a scaling solution) and D2 (a 
quintessence-dominated solution), each with wde > — 1, 
are stable. In these regions the initial conditions will 
determine the cosmological evolution. 

As we have seen, point F exists in two disconnected 
regions in the parameter space: a stable region where 
a < and an unstable region where a > 0. This can be 
understood as a fracture of the unstable regions in the 
parameter space, depicted in Fig. [2] In the standard 
quintom case, there are two overlapping regions corre- 
sponding to C2 and D2. When afj, — 2A > the points 
E and F are present and encroach upon this region. As 
one might expect, in both cases, when the exponent of 
the quintessence field is large, the scaling solution does 
exist (although it is not always stable). 



IV. COSMOLOGICAL IMPLICATIONS 

Having performed the complete phase-space analysis of 
the model, we can discuss the corresponding cosmologi- 
cal behaviour. A general remark is that the phenomenol- 
ogy is different from both the standard {a = 0) quintom 
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FIG. 1: The parameter space showing the stability of the 
critical points for (a) a — —3 and (b) a — 3, with 7 = 1. The 
light shaded area indicates the region where Dl is stable and 
the dark shaded area indicates the region where F is stable. 
There is a stable point with wde < — 1 at each region in the 
parameter space. The area marked in graph (b) indicates the 
region where C2 and D2 are stable. 



scenario [8j, |9[ and from the case of quintessence with 
mixed kinetic terms [l8j . In particular, there are late- 
time solutions describing both quintessence and phantom 
dominated universes, with the crucial observable quan- 
tity being the dark energy equation-of-state parameter 



« t i i i i i i i i i i i p 

1 2 3 4 5 6 

FIG. 2: Regions in the parameter space where the critical 
points exist but are unstable for a = 3 (with 7=1). 



iode (which is greater than —1 in the former case and 
less than —1 in the latter). 

The points A, B and E are unstable, and thus cannot 
be late-time solutions of the Universe. A corresponds to 
a matter-dominated universe, while B represents a dark 
energy-dominated one, with the stiff value ujde = 1- ^is 
a scaling solution. Finally, in one of the two disconnected 
areas in which F is physical, when ap, — 2A > and 
< II < 6, which corresponds to u>de > — 1, the point is 
always unstable. 

As usual, the stable critical points can be classified as 
field-dominated (Qde — 1) or scaling solutions (u>de = 
in the presence of matter). Since these features are in- 
compatible with observations, we deduce that the late- 
time behaviour of the scenario considered here cannot 
provide a satisfactory description of the Universe at the 
present time. However, the fact that dark energy has 
only recently begun to dominate the energy density, af- 
ter an extended matter-dominated era, suggests that the 
Universe may have yet to reach an attractor solution. In 
this case, one can consider our Universe as evolving to- 
wards a field dominant attractor, alth oug h the problem 
with initial conditions is still an issue |26[ . 

In order to treat the coincidence problem one must ex- 
plain why the present dark energy and matter density 
parameters are of the same order of magnitude. Alter- 
natively, one can ask why dark energy has only recently 
become dynamically important, which is a problem of 
initial conditions. This can be alleviated if pde scales 
with the dominant matter component, as the discrep- 
ancy between the initial energy densities of the various 
components of the Universe does not have to be so large. 

In contrast to the simple quintom scenario with a = 



7 



0, there is a stable scaling solution in this model, viz. 
point C2 7 which is an extension of the unstable point T 
in [§]. However, as the present model stands, there is 
no mechanism to effect a transition between the scaling 
regime and field-dominant solution without fine-tuning 
the initial conditions. 

Another novel feature of the extended quintom sce- 
nario is a stable field-dominated late-time solution with 
i«DE > — 1 (point D2) which is characterised by the fact 
that the contribution of the potential of the phantom field 
to the energy density is zero. Like C2, this is an extension 
of an unstable point in the simple quintom scenario. The 
deviation from iode = — 1 is 4 A 2 / [3(4 + a 2 )] so one re- 
quires cither the coupling to be large or the quintessence 
potential to be sufficiently flat for this solution to be 
compatible with observations. 

Let us make an important comment concerning C2 
and D2. As can be seen in Table HI both points have 
x ac > in the regions in which they are stable (we con- 
sider positive potential exponents X, fi > 0). However, 
since phantom fields generally do not roll down their po- 
tential slopes, one would normally expect that in a sta- 
ble late-time solution, the rate of change of the phantom 
field, that is x ac , should be negative (in a stationary solu- 
tion) or zero (in a static solution). Indeed, for the simple 
quintom case (a = 0) we do obtain x ac = 0. Interestingly 
enough, the positive coupling of the kinetic terms of the 
two fields makes the quintessence field drag the phantom 
field down its potential (i.e. the opposite to "normal" 
phantom behaviour). This possibility is a novel feature 
of the present scenario. 

The point Dl is the analogue of D2 and is stable for 
the larger part of the parameter space (it is always stable 
for a > 0, while for a < it can still be stable for > 
— a\/2). u>de depends only on the values of /x and a and 
approaches —1 if the potential is flat or the coupling is 
large, in a similar manner to the quintessence-dominated 
point, D2. Furthermore, it is interesting to note in this 
case there is a situation similar to the one described in the 
previous paragraph. In particular, for a < we observe 
that X(p c < 0, so, contrary to expectations, there is a 
stable late-time solution where the quintessence field rolls 
up its potential slope. Again, this is due to the coupling 
of the kinetic terms: the phantom field feeds the upward 
rolling of the quintessence field. 

From the discussion of the above two paragraphs we 
deduce that the kinetic coupling leads to new stable so- 
lutions. Additionally, it could lead to a very interesting 
evolution towards these late-time solutions. In particu- 
lar, if the kinetic terms are initially small, the effect of 
the coupling is downgraded, and the fields can move away 
from these stable points. With the increase of their ki- 
netic energies, the coupling becomes more important and 
they will start to move towards them. 

When point Dl is unstable, the critical point F, in one 
of the two disconnected areas for which it is physical, 
viz. 2/i + aX < 0, is stable. It corresponds to a universe 
completely dominated by dark energy, behaving like a 



phantom, with the EOS determined by the parameter H 
In this case, the quintessence and phantom fields combine 
to give a solution with wbe < — 1- This is the analogue 
of point F in [l8| and, like the unstable point E, does not 
exist in the standard quintom case. 

-0.6 | 1 1 1 1 1 1 1 , 1 

-0.7 - 




-1.2 T 



1 2 3 4 5 6 7 8 9 

N 

(a) 

0.2 i 1 1 1 1 1 1 1 1 1 




'0 1 2 3 4 5 6 7 8 9 

N 

(b) 

FIG. 3: Numerical integration of the system with (a) fi = 1, 
A = 1, a = 3 and (b) /i = 2.5, A = 3.5, a = 3 with initial 
conditions chosen so that the matter field (with wm — 0) 
initially dominates. A variety of plots with different initial 
conditions are shown: the colour-coding is for visual clarity. 
In case (a) points Dl and D2 are stable; one can see that there 
are solutions that cross the phantom divide to reach the late- 
time attractors. A similar situation is seen in case (b) where 
Dl and C2 are stable. 

In general, if a cosmological model possesses more than 
one stable critical point, amongst the possible late-time 
solutions it will reach the one in the basin of attraction 
of which it lies. Thus, the specific cosmological evolution 
depends heavily on the initial conditions. In our scenario, 
there are regions of the parameter space where there is 
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both a stable quintessence-like and a stable phantom- 
like critical point. Therefore, the EOS can always re- 
main in the quintessence or the phantom regime, cross 
the phantom divide from below to above, or cross the 
phantom divide from above to below, with the last possi- 
bility mildly favoured by observations. It must be noted, 
however, that in the absence of constraints on the cou- 
pling a and the potential exponents, a late-time solution 
with wde < — 1 seems likely as there is a stable phantom- 
like critical point in every region of the parameter space. 
This fact is an advantage of the model. 

As we know, the standard quintom model (8l-fl0j ex- 
hibits only one (phantom-like) stable late-time solution, 
and in order to acquire the simultaneous existence of 
one quintessence and one phantom stable late-time so- 
lutions, one has to extend to significantly generalized 
quintom models [24[. On the other hand, the model at 
hand presents this behavior more economically. Further- 
more, note that the extended scenario of this work ex- 
hibits a novel situation, viz. the presence of more than 
one stable attractor (points C2 and D2) corresponding to 
quintessence-like solutions. Thus, according to the spe- 
cific parameter values, and for suitable initial conditions, 
the late-time behaviour of the Universe can be described 
by one of the quintessence-like states, with D2 possessing 
a more physical EOS compared to C2 (see Table HT| . 

In order to present the late-time behaviour transpar- 
ently, in Fig. [3] we depict the results of a numerical in- 
tegration for two choices of parameter values. In graph 
(a) the phantom stable point is Dl and the quintessence 
one is D2, while in graph (b) the phantom stable point 
is also Dl but the quintessence one is C2. If the system 
avoids the stable phantom solution, it will tend toward 
one of the two quintessence-like solutions, depending on 
the values of the parameters, a, [i and A. 

We have to mention that for a particular subclass of 
evolutionary trajectories, depending on the parameter 
values and initial conditions, the dark energy EOS be- 
comes infinite at a finite time t*. However, this diver- 
gence in wde is not accompanied by a divergence in the 
scale factor, its time-derivative or in the dark energy den- 
sity and pressure, but is an expected result arising from 
the nullification of pde, possible in all phantom models. 
Thus, technically, it does not correspond to the Big Rip 
discussed in the literature 12711 . but rather to the new sin- 
gularity family discussed in [28j . One could consider only 
those models in which the total energy density remains 
positive (where the potential is sufficiently positive or the 
phantom kinetic energy not too large), to avoid the cor- 
responding energy condition violation. Finally, note that 
while a negative coupling value a facilitates the nullifica- 
tion of pde, a positive a obstructs it. 

As one can observe in Fig. [3l the dark energy EOS can 
change significantly before the system finally reaches the 
critical points. It is therefore of great interest to investi- 
gate whether the situation described above is compatible 
with the stringent bounds provided by cosmological ob- 
servations. To this end, in Fig. [3] we replot the evolution 



-0.8 



-0.85 - 




0.5 1 1.5 2 2.5 3 3.5 4 

Z+1 



(b) 

FIG. 4: The evolution of the dark energy EOS with respect 
to redshift [z + 1 = a(to) /o(t)]. Parameter values, initial con- 
ditions and colour-coding are as in Fig. [4] z — is chosen to 
correspond to the point at which Q.de = 0.728. Shaded re- 
gions are excluded by Wilkinson Microwave Anisotropy Probe 
7 yr data+BAO+ffo+SN at 68% confidence level Q. Using 
the parameterisation iude = Wo + w a z/(z + 1), the limits are 
too = 0.93 ± 0.13 and w a = -0.41±g;^. 



of Wr>E of Fig. [3l using the redshift z [z + 1 = a(t )/a(t)] 
as the independent variable, presenting also the obser- 
vationally excluded regions from [29j . up to z = 3. In 
the upper graph, both potentials are relatively fiat and 
the values of wde at the critical points Dl and D2 are 
close to —1. This means that, despite the range of be- 
haviour exhibited for different initial conditions, the EOS 
does not deviate much from —1 when the trajectories are 
approaching the attractors at recent times. These fea- 
tures make the evolutions of Fig. |4ja) compatible with 
observations. 

The case in Fig. [5£b) is quite different. As mentioned 
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previously, it can be seen that the scaling point C2 is 
ruled out as a late-time solution. Also, for both C2 and 
Dl to be stable, A needs to be larger, leading to a more 
negative EOS for Dl. As the trajectories start to con- 
verge on this value, wue evolves more quickly than Fig. 
HJa), leading to an EOS that is disfavoured by observa- 
tions at z w 0. 

The behaviour of Fig. 0] is fairly typical. In the case 
where the dark energy contribution to the energy content 
of the Universe is initially negligible, the value £Ide ~ 
0.7 is obtained only when the solutions start to converge 
to the critical points. Thus, if the value of wde at the 
critical points is closer to —1, the variation in the EOS 
is smaller at recent times, and closer to the observational 
values. 

In order to perform a complete comparison with the 
data, it is necessary to take into account the effect of 
dark energy perturbations. In particular, the late-time 
integrated Sachs- Wolfe effect, which affects the CMB on 
large angular scales, is sensitive to the presence of dark 
energy perturbations [22[ • Perturbations in the standard 
quintom model have been investigated in [30j . where it 
was found that the parameter space that allowed the EOS 
to cross the phantom divide was enlarged when the per- 
turbations were included. As the present scenario ex- 
hibits behaviour similar to simple quintessence and quin- 
tom models, we expect that many features of these sce- 
narios will be retained when the perturbations are taken 
into account. However, as we have seen, an important 
difference between the coupled and uncoupled cases is 
the presence of a source term proportional to a that can 
drag the quintessence (phantom) field up (down) its po- 
tential. This would introduce extra terms into the pertur- 
bation equations which could lead to observable effects. 
Another interesting possibility is that the behaviour of 
isocurvature perturbations will be affected by the trans- 
fer of energy between the two fields. The full analysis 
of the perturbations in the model at hand is beyond the 
scope of the current investigation and will be treated in 
a future work (3lj . 



cal points E and F that are not present in the simple 
phantom model, the latter of which is analogous to that 
responsible for the phenomenon of assisted quintessence 
in the case with two canonical fields. Here, however, 
the combined effect of the fields is to give phantom-like 
behaviour, with wde < — 1. As well as this, there is an 
extension of the phantom-dominated asymptotic solution 
in the standard quintom model, which is stable for a large 
region of the parameter space. In all, there is a stable 
solution with u>de < —1 and £Ide = 1 in every part 
of the parameter space, even when the quintessence-like 
points are stable. We have discussed the cosmological 
implications of these results in Sec. HVl 

A characteristic feature of the presence of the kinetic 
coupling is the presence of stable solutions in which the 
quintessence field can be dragged up its potential, in the 
opposite direction to the standard case. This occurs 
when the coupling is negative, where the quintessence 
equation of motion acquires a positive source term. Sim- 
ilarly, for a positive coupling the phantom-field equation 
of motion obtains a negative source term, which causes 
the field to roll down its potential. This could affect the 
evolution of the dark energy perturbations and lead to 
observational differences between the kinetically coupled 
quintom scenario and other models, offering a relatively 
safe signature. 

The use of canonical and phantom fields to drive the 
acceleration of the Universe in phenomenological mod- 
els is not necessarily representative of the fundamental 
mechanism behind dark energy, in which non-trivial fea- 
tures such as non-canonical kinetic terms and kinetic cou- 
plings may be important. The scenario of the present 
work is interesting as the kinetic coupling gives rise 
to behaviour different from that of the standard quin- 
tom model. In particular, both quintessence-like and 
phantom-like late-time solutions, as well as solutions that 
cross the phantom divide iude = — 1 during the evolution, 
are possible. Thus, such scenarios could be a candidate 
for the description of dark energy. 



V. CONCLUSIONS 

We have considered an extension of the quintom model 
of dark energy in which the kinetic terms of the phan- 
tom and the canonical scalar field are coupled by a term 
in the Lagrangian of the form ^g^d^d^a. In order to 
study the asymptotic behaviour of the model and to fa- 
cilitate comparison both with the standard quintom sce- 
nario and the two-field quintessence model with mixed 
kinetic terms, we have performed a phase-space analysis 
of the corresponding dynamical system. 

We find that the kinetic interaction allows for the pos- 
sibility of stable critical points similar to those found 
in quintessence scenarios, including field-dominated so- 
lutions with u>de > — 1 and solutions displaying scal- 
ing behaviour. Additionally, there exist two new criti- 
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Appendix: The linearised perturbation matrix 

The 4x4 matrix Q of the linearised perturbation equa- 
tions of the system (fTS)) reads: 

Qn Q12 Q13 Qi4 

q _ Q21 Q22 Q23 Q24 

Q31 Q32 Q33 Q34 

_ Qu Q42 Q43 Q44 . 
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with 



hi = -3 + |a;0 C (2 - j)(2x^, c + ax ac ) + T c , 
?12 = §2^ c (2 - 7)(-2x (TC + ax^c), 



Ql3 — Vtj>c 



hi = Vac 



(4 + a 2 ) 
(4 + a 2 ) 



SjXrjtc 
SjX^c 



Q21 
Q22 

Q23 
Q24 



Q31 — y<pc 



\x ac (2 - r )){2x (f)C + ax ac ), 

-3 + \x ac (2 - "f)(~2x ac + ax^ c ) + T c 



2aV6\ 
(4 + a 2 ) 

(4 + a 2 ) 



37X 



37X 



[(2 - 7)(2x^ c + ax ac ) 



Q32 = \y<$>c{2 - 7)(-2x ac + ax jc), 



Q42 
Q43 
Q44 



where 



Q33 = -y 2^ x 4>c - 372/^c + T c 
Q34 = S^cVac, 



\v<yc{2 - j)(2x 4>c + ax ac ), 
yac |(2 - 7)(-2x ac + ax<p c ) - 



2 l^*E<rc 



^IVac + T c 



2x, - 2x ac + 2ax cf)C x crc + 



(A.l) 



+7(1 - x\c + %lc - uXtfycXvc - y\ c - y 2 ac ) . (A.2) 
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